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ABSTRACT 

Galaxies exhibit a sequence of various morphological types, i.e., the Hubble sequence, and they are ba- 
sically composed of spheroidal components (elliptical galaxies and bulges in spiral galaxies) and disks. It 
is known that spheroidal components are found only in relatively massive galaxies with M — lO^^^^^A/©, 
and all stellar populations in them are very old, but there is no clear explanation for these facts. Here we 
present a speculative scenario for the origin of the Hubble sequence, in which magnetic fields ubiquitously 
seen in galaxies have played a crucial role. We first start from a strange observational fact that magnetic 
field strengths observed in spiral galaxies sharply concentrate at a few microgauss, for a wide range of 
galaxy luminosity and types. We then argue that this fact and the observed correlation between star 
formation activity and magnetic field strength in spiral galaxies suggest that spheroidal galaxies have 
formed by starbursts induced by strong magnetic fields. Then we show that this idea naturally leads to 
the formation of spheroidal systems only in massive and high-redshift objects in hierarchically clustering 
universe, giving a simple explanation for various observations. 

Subject headings: cosmology: theory — galaxies: formation — galaxies: magnetic fields — stars; 
formation 



1. INTRODUCTION 

It is generally believed that the hierarchical structure 
today seen in the universe is generated by gravitational 
instability of cold dark matter (CDM) whose density is 
about 10 times higher than that of baryonic matter, and 
galaxies form in virialized dark matter haloes by cooling of 
baryonic gas and subsequent star formation. The charac- 
teristic mass range of galaxies is about \(f'~^'^MQ, which 
can be understood as a mass range in which baryonic gas 
can cool sufficiently within a dark matter halo to form a 
galaxy (e.g., Blumcnthal et al. 1984 for a review). How- 
ever, spheroidal components reside only in relatively mas- 
sive galaxies with M = 1{)^^~^'^Mq and less massive galax- 
ies are mostly irregular types or late-type spiral galaxies 
(Burstein et al. 1997; Hunter 1997). [There is a pop- 
ulation called as dwarf spheroidal galaxies, but they are 
generally considered to be different from giant elliptical 
galaxies and bulges. They are more similar to dwarf ir- 
regulars (Bender et al. 1992; Mao & Mo 1998).] Ellipti- 
cal galaxies and bulges appear to form a very uniform, 
old stellar population with very little scatter in metal- 
licity and star formation history, and they are generally 



considered to have formed by intensive starbursts at high 
redshifts (e.g., Renzini 1999). The trigger mechanism for 
starbursts is unknown, and it is often assumed to be merg- 
ers of disk galaxies when galaxy formation is modeled 
in the context of hierarchical structure formation in the 
CDM universe (KauflFman, White, & Guiderdoni 1993; 
Baugh, Cole, & Frenk 1996). However, there is no ev- 
idence for spheroidal galaxy formation through mergers, 
and the merger hypothesis does not give a clear explana- 
tion for the fact that all spheroidal galaxies are very old 
and only in massive galaxies. These trends appear to run 
opposite to the expectation of structure formation in the 
CDM universe, in which smaller objects form earlier and 
then massive objects later. This is currently considered as 
a major challenge for galaxy formation in the CDM uni- 
verse (e.g., Renzini 1999), and it is interesting to seek for a 
mechanism which triggers starbursts only in massive and 
high-redshift objects. In this letter we show that a possible 
strong dependence of star formation activity on magnetic 
field strengths, which is suggested by observations of mag- 
netic fields in local galaxies, gives a candidate for such a 
mechanism. 
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2. GALACTIC MAGNETIC FIELDS AND STAR 
FORMATION ACTIVITY 

Interstellar magnetic fields in galaxies are strong enough 
to significantly affect interstellar gas dynamics on both 
global scales characteristic of galactic structure and small 
scales characteristic of star formation (see, e.g., Vallee 
1997; Zweibcl & Heiles 1997 for a review). Fitt and 
Alexander (1993; hereafter FA) derived strengths of 
volume-averaged magnetic fields for 146 spiral galaxies, 
which is a complete sample including various types from 
Sa to irregular galaxies, and they found a strange result: 
the dispersion in strength is surprisingly narrow within a 
factor of less than 2 with an average strength of 3 /uG, in 
spite of a wide range of absolute radio luminosity of galax- 
ies (almost 4 orders of magnitude). Why is magnetic field 
strength so uniform at ^ 3/iG for various types of galax- 
ies with such a wide range of luminosity? Clearly there 
should be a physical reason for this strange fact. 

Vallee (1994) investigated the relation between the mag- 
netic field strength and star formation activity by us- 
ing 48 galaxies out of the FA sample, and found that 
there is a correlation between magnetic field strength (B) 
and star formation efficiency (SFE) as B oc (SFE)* with 
i = 0.13 ± 0.04. [The same relation was found also for star 
formation rate (SFR) and B.] From this result, one may 
conclude that star formation activity docs not strongly af- 
fect the strength of magnetic fields. However, it is also 
likely that star formation rate is affected by magnetic 
fields. If we take the result of Vallee in this context, the 
observed correlation suggests that star formation activ- 
ity depends quite sensitively on B as SFE oc B"^ with 
q = 1/i ~ 7.7 ± 2.6. Here we point out that, if this 
is the case, the surprisingly narrow dispersion in galactic 
magnetic fields is naturally explained by an observational 
selection effect. If the strength of magnetic field in an 
ob ject is significantly smaller than a few /iG, the star for- 
mation activity becomes much lower than in typical spiral 
galaxies, and such object cannot produce enough stars to 
be observed as a galaxy. On the other hand, if magnetic 
fields are stronger than a few /iC, such objects would ex- 
perience strong starburst and all interstellar gas will be 
converted into stars. If the star formation time scale is 
shorter than the dynamical time scale for disk formation, 
they will become present-day spheroidal systems. Mag- 
netic fields in elliptical galaxies 'atc difficult to measure 
since there are no relativistic electrons to illuminate the 
magnetic fields by synchrotron radiation. Then the rea- 
son why the observed dispersion of magnetic fields is quite 



small can be understood by an observational selection ef- 
fect; we cannot observe or measure magnetic field strength 
when the strength is significantly different from the typi- 
cally observed value of ~ 3/xG. 

Therefore the observed correlation between B and SFE, 
as well as the surprisingly narrow dispersion of measured 
values of magnetic field strength, well suggests a strong 
dependence of star formation activity on magnetic field 
strength. This idea is further supported by an amorphous 
spiral galaxy M82, which has abnormally large magnetic 
field of ~ 10/LtG in the sample of FA. M82 is actually 
known as an archetypal starburst galaxy, and the strong 
field strength suggested by FA is caused by the very strong 
field (^ 50y^G) in the nuclear region of this galaxy, where 
intense star formation is underway (Klein, Wielcbinski, & 
Morsi 1988; Vallee 1995a). Unfortunately our knowledge 
for physics of star formation is poor and it is not clear from 
a theoretical point of view why SFE depends so strongly 
on B, but it is not unreasonable. In fact, at least one effect 
is known by which magnetic fields help star formation: the 
magnetic braking (e.g., McKee et al. 1993 for a review). 
In order for a protostellar gas cloud to collapse into stars, 
a significant amount of initial angular momentum must 
be transported outward. Magnetic fields play an impor- 
tant role for this angular momentum loss by Alfven waves 
launched into ambient medium. Strongly magnetized ob- 
jects with large scale high gas density would be followed 
by magnetic braking of individual regions with angular 
momcntiim losses, ending in outgoing Alfven waves, gas 
collapse, and star formation on a large scale. 

It should be noted that the above hypothesis does not 
contradict with the well-known Schmidt law of star for- 
mation. Vallee (1995b, 1997) found a relation between 
magnetic field strengths and interstellar gas density on a 
scale of a whole galaxy as _B oc n'^ with k = 0.17 ± 0.03. 
By using this relation and S-SFR relation, we find SFR 
oc n'^/' with k/i = 1.3 ± 0.1. This relation agrees well 
(Vallee 1997) with some direct estimates of the relation 
between SFR and n (e.g., Kennicutt 1989; Shore & Ferrini 
1995; Niklas & Beck 1997). At the same time, these rela- 
tions suggest that the thermal energy of interstellar matter 
is not in equipartition with the magnetic field energy or 
cosmic ray energy (Kamaya 1996). It is possible that what 
is physically affected by magnetic fields is gas density, and 
then the gas density determine SFE by the Schmidt law, 
producing the observed S-SFE relation. 

Anyway, the correlation between B and SFE has ac- 
tually been observed, and it is natural to consider this 
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relation holds also in high redshift objects. [The chance 
probability to observe such a correlation from an uncorre- 
lated parent population is only 12% (Vallee 1994).] There- 
fore the following analysis is valid unless some unknown 
processes violate this empirical relation in high redshift 
objects. 

3. MAGNETIC FIELDS IN HIERARCHICAL 
STRUCTURE FORMATION 

As discussed above, it is suggested that spheroidal com- 
ponents of galaxies formed by starbursts induced by strong 
magnetic fields. In the following we consider the magnetic 
field generation and formation of spheroidal systems in 
the framework of the standard cosmological structure for- 
mation in the CDM universe, and show that spheroidal 
systems form only in relatively massive objects at high 
redshifts. Recently it is discussed that galactic magnetic 
fields are generated at the stage of collapse of proto- 
galactic clouds (Kulsrud et al. 1997). Although there 
are some other scenarios as the origin of galactic mag- 
netic fields, such as earlier field generation on cosmologi- 
cal scales or dynamo amplification after galaxies form (see, 
e.g., Zweibel & Heiles 1997), we assume here that mag- 
netic fields are generated during collapse of protogalactic 
clouds. When a dark matter halo decouples from the ex- 
pansion of the universe and collapses into a virialized ob- 
ject, the gravitational energy of baryonic gas is converted 
into turbulent motions and thermal energy by generated 
shocks. It can be shown that a very weak magnetic field 
(~ 10~^^ G), which has been generated by thermoelec- 
tric current before collapse (i.e., the battery mechanism), 
is amplified up to a strength nearly in equipartition with 
turbulent energy (Kulsrud & Anderson 1992; Kulsrud et 
al. 1997). The time scale for equipartition is given by 

r/u, where r is the length scale of the system and v 
the turbulent velocity. (The essence of this result can be 
understood by a dimensional analysis. The equation of 
magnetic field generation is dB/dt = rot(v x B), where v 
is the velocity field of fluid, and when we estimate rot by 
the inverse of the system scale, it is clear that the field 
evolution time scale is given hy r/v.) If we take the 
radius and three-dimensional velocity dispersion of a viri- 
alized halo as r and v, it is straightforward to see that 
this time scale is given by the dynamical time scale of the 
halo. Then we can estimate the strength of magnetic fields 
by the turbulent energy density of baryonic gas in a col- 
lapsed object, which is determined by using the well-known 
spherical collapse model (e.g., Peebles 1980). For a dark 



halo with mass Mh and formation redshift z, the turbu- 
lent energy density becomes e_B ^ {SilBMhv"^ /Sir^or^) = 

6.4 X 10~13/j8/3^j2/3^^ _|_ ^)4^^ gj.g ^Y^^j.^ jg i-jjg 

Hubble constant normalized at 100 km/s/Mpc, M12 = 
Mh/{10^^MQ), and and fig are the baryon and mat- 
ter density in units of the critical density in the uni- 
verse. In the second expression above, we have assumed 
the Einstein-de Sitter universe (fio — 1), but extension 
to other cosmological models is easy. According to the 
theory of Kulsrud and Anderson (1992), we assume that 
one-sixth of the turbulent energy is converted into mag- 
netic field energy and then the equipartition magnetic field 
becomes B - IM^^^mI/,^ {1 + zffl]^^ fiG. The observa- 
tional properties of disk galaxies including our Galaxy are 
well understood if they are considered to have formed at 
z ~ 0-1 (Mao & Mo 1998), and hence this theory of mag- 
netic field generation gives a roughly correct strength for 
our Galaxy with Mh IO^^Mq. 

4. EMERGENCE OF THE HUBBLE SEQUENCE 

Then spheroidal galaxies are expected to form at high 
redshifts in massive dark matter halos, because magnetic 
field becomes stronger with increasing mass and redshift as 
B cx il/^^'^(l + z)^. In order to discuss more quantitatively, 
we equate the previously defined star formation efhciency 
(in §2) tov= {tsF)~'^, where tgF is the time scale on which 
interstellar gas is converted into stars. The star formation 
rate M* in a galaxy is given by M* = Mgas/isF, where 
Mgas is the mass of interstellar gas in the galaxy. As men- 
tioned earlier, the observed relation between B and v is 
v (X B"^ with q ^ 7.7±2.6. We set the normalization of this 
relation by SFE at our Galaxy. The star formation rate in 
our Galaxy is about a few Mq / yr and mass of interstellar 
gas in the disk is about 6 x IO^Mq at present, suggesting 
that tsF ~ 2 Gyr at B ^ 3/iG (e.g., Binney & Tremaine 
1987). We also normalize the strength of S to 3 /j,G for a 
halo whose baryonic mass is 6 x IO^^Mq and whose viri- 
alization occurred at z = 1, supposing our Galaxy. Then 
we can calculate the ratio of baryonic dynamical time to 
tsF, r = idyn/isF, which Can be considered as a criterion 
for spheroidal formation. In calculation of idynj we use a 
baryon density A~"^ times higher than the value at viri- 
alization, considering the contraction of baryons due to 
cooling and dissipation, where A ~ 0.05 is typical dimen- 
sionless angular momentum of dark haloes (Warren et al. 
1992). If r ^ 1, spheroidal systems are expected to form. 

In Figure |l|, we have plotted the contour of F = 1, 10, 
and 100 by thin solid lines in a plane of baryon mass of a 
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collapsed object (Mbaryon) and its formation redshift, as- 
suming q = 5. (Observed mass of galaxies is considered to 
be mainly baryonic mass.) For reference, we have also plot- 
ted the contour of magnetic field strengths in Fig. ||. We 
have used a flat, A-dominated cosmological model, where 
A is the cosmological constant, with the standard cos- 
mological parameters of (/i, flo, ^^A, cs) = (0.7, 0.3, 0.7, 1). 
This universe is well consistent with various observations 
such as ages of globular clusters, high redshift type la su- 
pernovae (Perlmutter et al. 1998), cosmic microwave back- 
ground (Bunn & White 1997), and abundance of clusters 
of galaxies (Kitayama & Suto 1997). In the following we 
discuss the emergence of the Hubble sequence using Fig. 
1^. Objects lying above the thin solid lines are expected 
to form spheroidals, but we have to check how such ob- 
jects are typical in the context of cosmological structure 
formation in the CDM universe. For this purpose we have 
plotted mass scale of n-a fluctuation as defined in Blumen- 
thal et al. (1984), by dashed lines. The lines are defined by 
na{AIh, z) = Sc, where a{Mh,z) is the root-mean-square 
of density fluctuation predicted by linear theory in the 
CDM universe at scale of Mh and at redshift z (Peacock 
& Dodds 1994; Sugiyama 1995), and Sc ~ 1.69 is a criti- 
cal density contrast at which an object virializes (Peebles 
1980). The three dashed lines correspond to n = 1, 2, 
and 3, and these lines represent typical mass scales of col- 
lapsed objects as a function of z. (For larger n, the cosmo- 
logical abundance of objects is statistically suppressed as 
oc e~" /^.) In the region where the dashed lines are above 
the solid lines, spheroidal galaxies can form as cosmolog- 
ically typical objects. Figure |l| shows that objects with 
A^baryon > 10^^^^ Mq with ^ 2-3(7 fluctuation will form 
spheroidal systems because F 1, at high redshifts of z 
3-10 with B 50^G (see also Fig. ||). It is interesting to 
note that this magnetic field strength is roughly the same 
with that in the starbursting nucleus of M82, as mentioned 
in §2. Spheroidal formation with Afbaryon < 1O^~^"M0 
is inhibited because of low F for cosmologically typical 
objects. On the other hand, if Afbaryon is larger than 
~ 1O^^M0, formation of galaxies is inhibited by too long 
cooling time (icooi) compared to the Hubble time at each 
redshift {tn), as shown by dotted lines which are contours 
oficooi/^H (Blumenthal et al. 1984). Therefore, our model 
can explain the observed mass range {10^^~^'^Mq) and old 
stellar populations seen in spheroidal galaxies. The region 
for spheroidal galaxy formation is schematically indicated 
in Figs. |l| and || by the thick solid line. After spheroidal 
galaxy formation, gas accretion onto some of them is pos- 



sible at lower redshifts. Since such gas accreting recently 
is located below the thin solid lines in Figure 0, it results 
in disk formation of spiral galaxies at z < 1. On the other 
hand, if recent gas accretion is negligible, such object will 
be seen as present-day elliptical galaxies. 

Since solid lines (constant F) and the dashed lines {n-a 
lines) are approximately parallel in the region of spheroidal 
formation, our scenario for the origin of the Hubble se- 
quence predicts that elliptical galaxies should lie on a con- 
stant n-a density fiuctuation with n ~ 2-3, and with de- 
creasing n, the type of galaxies becomes later in the Hub- 
ble sequence, i.e., from early to late spiral, and then into 
irregular galaxies at ^ Icr fiuctuation. In fact, this trend 
is exactly what has been observed in galaxies (Blumen- 
thal et al. 1984; Burstein et al. 1997), giving a further 
support for our scenario. Most properties or relations ob- 
served in present-day galaxies, such as the Tully-Fisher 
or Faber- Jackson relation, existence of the fundamental 
plane for galaxies and their distribution on it, can be ex- 
plained by formation of early type galaxies from higher 
n-a fiuctuations than later Hubble types (Burstein et al. 
1997). The density- morphology relation, which is the cor- 
relation between galaxy types and number density, is also 
explained; higher n-a fluctuations occur preferentially in 
denser regions destined to become rich clusters, and hence 
one expects to find more ellipticals there, as is observed 
(Blumenthal et al. 1984). Higher n-cr objects are expected 
to show stronger spatial clustering than lower ones (e.g.. 
Mo & White 1996), and it is consistent with the stronger 
clustering observed for elliptical galaxies than late type 
galaxies (Loveday et al. 1995). Various observations sug- 
gest that giant galaxies formed at higher redshifts of z > 3, 
and were then followed by a sequence of less and less mas- 
sive galaxies forming at lower and lower redshifts, lead- 
ing down to the formation of dwarfs at recent (z < 0.5) 
(Fukugita, Hogan & Peebles 1996; Sawicki, Lin, & Yee 
1997). The proposed scenario gives an explanation for this 
trend which is sometimes termed as "downsizing" , oth- 
erwise it seems opposite to the expectation in the CDM 
universe. 

5. DISCUSSION & CONCLUSIONS 

For very low mass objects, star formation is strongly 
suppressed by the absence of physical triggers when mag- 
netic fields are weak, and this may give an explanation for 
the fact that the faint end of luminosity function of galax- 
ies is much fiatter than expected from the mass function of 
dark haloes (Kauffmann et al. 1993; Baugh et al. 1996). 
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Because of the strong dependence of star formation activ- 
ity on magnetic fields, stellar luminosity of galaxies would 
quite rapidly decrease with decreasing mass of galaxies. 
Then it is expected that, in the faint end of galaxy lumi- 
nosity function, the mass of galaxies does not so change 
compared to the change in luminosity. Since the number 
density of objects is determined by mass of objects, the 
faint end is expected to be relatively flat, as observed. Ir- 
regular galaxies or late-type galaxies are expected to show 
a wide range of star formation activity within a narrow 
range of galaxy masses depending on the magnetic field 
strength in them, and in fact such a trend has been ob- 
served (Hunter 1997). 

It is well known that galaxy interactions or mergers in- 
duce intensive starbursts, but the mechanism which trig- 
gers starbursts in galaxy interactions is still poorly known. 
Mergers are followed by cloud collisions leading to high gas 
density, strong turbulences and strong magnetic fields, and 
ending in gas collapse and star formation in clouds. There- 
fore the hypothesis presented in this letter, i.e., strong 
dependence of star formation activity on magnetic field 
strengths, may also be important in the merger-induced 
starbursts. Some fraction of elliptical galaxies may have 
formed by such a process, but we have argued that the 



observed trends of ellipticals, i.e., being massive and old, 
originate mainly from properties of gravitationally bound 
objects in the standard theory of cosmological structure 
formation. 

We have presented a new idea for the origin of the Hub- 
ble sequence of galaxies. The key of this scenario is the 
strong dependence of star formation activity on average 
magnetic fields in galaxies, which is just a speculation from 
a theoretical point of view, but motivated by several ob- 
servational facts about magnetic fields in nearby galaxies. 
This only one speculation provides a simple explanation 
for the surprisingly narrow dispersion in the magnetic field 
strengths observed in spiral galaxies, and also for most 
properties of galaxies seen along the Hubble sequence. To 
verify this "magnetic galaxy formation" scenario as the 
origin of the Hubble sequence, it is indispensable to con- 
firm observationally that starbursts are actually triggered 
by stronger magnetic fields. Comprehensive measurements 
of galactic magnetic fields in larger samples of galaxies are 
necessary, especially for starburst galaxies at local as well 
as at high redshifts. 

The author would like to thank T. Kitayama, M. 
Shimizu, and an anonymous referee for useful comments 
and discussions. 
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Redshift z 

Fig. 1. — Thin solid lines show contours of the ratio of dynamical time scale to star formation time scale, tdyn/tsF = 1, 
10, and 100, as functions of baryon mass and formation redshift of cosmological objects. When tdyn/isF ^ 1, spheroidal 
systems are expected to form. Dashed lines show the mass-rcdshift relation for objects with a given n-a fluctuation for 
n = 1, 2, and 3, where a is the root- mean-square of the density fluctuations in the universe at a given mass scale. The 
value of n represents statistical abundance of collapsed objects. Dotted lines are contours of the ratio of cooling time to 
the Hubble time at given redshifts, for tcooi/ii? = 1) 3.2, 10. Galaxy formation is inhibited for objects with tcooi/tH > 1- 
The region in which spheroidal galaxies form is schematically shown by the thick solid line. A A-dominated flat universe 
with cosmological parameters of {h, fio, ^^a, f s) = (0.7, 0.3, 0.7, 1) is assumed, where A is the cosmological constant. 
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Redshift z 

Fig. 2. — The same as Fig. |l|, but for the contour of magnetic field strengths in gravitationally bound objects. The 
indicated values associated to thin solid contours are magnetic field strengths in /iG, and dotted lines are contours with a 
step of 0.5 in logi?. The magnetic field strength is normalized to 3 at Mbaryon = 6 x 10^^ Mq and 2 = 1 supposing our 
Galaxy, and this normalization is very close to the equipartition field with the turbulent energy density of protogalactic 
clouds (see text). The region in which spheroidal galaxies form is schematically shown by the thick solid line (see Fig. 



